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ABSTRACT 



O-t 

I ■ Context. New instrumental capabilities and the wealth of astrophysical information extractable from the near-infrared wavelength 

O ' region have led to a growing interest in the field of high resolution spectroscopy at 1-5 yum. 

^ , Aims. We aim to provide a library of observed high-resolution and high signal-to-noise-ratio near-infrared spectra of stars of various 

C/3 . types throughout the Hertzsprung-Russell diagram. This is needed for the exploration of spectral features in this wavelength range 

■ and for comparison of reference targets with observations and models. 
Methods. High quality spectra were obtained using the CRIRES near-infrared spectrograph at ESO's VLT covering the range from 
0.97 yum to 5.3 yum at high spectral resolution. Accurate wavelength calibration and correction for of telluric lines were performed by 

" , fitting synthetic transmission spectra for the Earth's atmosphere to each spectrum individually. 

• Results. We describe the observational strategy and the current status and content of the library which includes 13 objects. The first 

■ examples of finally reduced spectra are presented. This publication will serve as a reference paper to introduce the library to the 
I community and explore the extensive amount of material. 

I Key words. Atlases - Infrared: stars - Stars: atmospheres - Line: identification - Methods: data analysis - Techniques: spectroscopic 

(N 

o: 

' 1 . Introduction tion in the NIR compared to the visual range allow us to access 

■ more distant stars. 
After a long-term reliance on the visual spectral region, it has r., ^ , , , ■ r ■ • . • . 
, 1 1 J » » J • J J *u To rully understand the information contained in the 
been clearly demonstrated in the past two decades that the near- ■' , ■ ■ ■ 

, ., ■ i , , , 11 1 1 J- 1 • NIR spectral range, high-resoIution spectroscopy is manda- 

and mid-infrared spectral ranges will play a leading role in many ^ • i , , • • i , i , 

.... c !^ 1 u c • J J tory. Considerable progress in instrument development has been 

prune areas of astronomical research for the coming decades. , . , , , • • , 

rrt ' Sn,- . f . ui u c ■ c achieved over the past twenty years, resulting in the current 

^« . This part of the spectrum enables observation of a universe of r>z ■ i ^- i xTt^. rr.i ■ i 

" " „ if, ■ • 1 J- J- 1 1 . . J J availability of elncient spectrographs for the NIR. The lack of 

cool phenomena, including disks, planets, or the extended ... , , • , i • i • , 

c 1 J . T-i . • . V c ^ hiRh quality (high spectral resolution and high signal-to-noise 

atmospheres of evolved stars. Electronic transitions of atoms, • r, Ti • i t^. i i 

, . ^, r 1 1 ^ ^- 1. ^- ^ -^^ ratio - S/N) stellar spectra in the IR has meant that spectro- 

and in the case of molecules rotation-vibration transitions, pro- . , , i • i ^ 

, ,. • • r J /KTTn\ ■ J- . £ scopists have concentrated their resources on the strongest fea- 

duce lines in the near-infrared (NIR), among them indicators for ^ , , , , , • • i t^. • i 

. I f A tures observed at lower spectral resolution in the IR or in other 

obtaining isotopic ratios and resonance lines of .s-process and , i • . , ,• i , , 

? , ^ ^ ^, 1 , , ,1 , wavelength regions. As a result, line data for many elements, but 

rare earth elements. For the coolest stars, the optical spectrum • ■ • , , , i • , 

. , • , J , 1 1 1- 1- . . . u- u especially the post-iron group elements, are severely lacking and 

is donunated by molecular lines, limiting the extent to which , , ■ • , ■ , 

, . 1 u J . • J • . ui J- A* those that exist are predominantly the result of theoretical calcu- 

chemical compositions can be determined owing to blending. At , . , , , /„r n i ^^^^ 

■ c J /Tr.\ 1 . 1 ■ J : lations. The atlases of the solar spectrum (Wallace et al. 1996, 

infrared (IR) wavelengths there are spectral windows that are ^ .,, , , , ,/ . . , , , 

, . , \, f 11 , ^- u V c u Delbouille et al. 1973) and the K giant Arcturus (Hinkle et al. 

relatively free of molecular absorption, where lines from heavy , i , .c n .c i - i 

1 / 1 . J r- J A cc ^ c ^- 1995) illustrated the value of a stellar reference spectrum at high 

elements are located. Furthermore, the reduced effects of extinc- , , . i, - i o/xt 

spectral resolution and high S/N. 

At high resolving power {R - 100 000), an inventory of com- 



><: 



* Based on observations collected at the European Organisation for , _ „_ „ , ,t i-. h i- 

Astronomical Research in the Southern Hemisphere, Chile (084.D- Pl^te NIR stellar spectra across the Hertzsprung-Russell diagram 

0912, 085.D-0161, 086.D-0066, and 087.D-0195. The spectra pre- (HRD) is presently unavailable. Spectra of both high spectral 

sented in Figs. 3 to 15 are available at the CDS via anonymous ftp to resolution and high S/N are key demands for a list of primary 

cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi- science cases in addition to the above-mentioned exploration of 

bin/qcat?J/A+A/ weak atomic lines of hardly studied elements. In the case of 
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molecular lines, both prerequisites are needed to study isotopic 
ratios, which open a window to the nuclear fusion processes in 
the stellar intererior. Another topic that can be investigated with 
the help of such data are small velocity shifts between high and 
low excitation lines or asymmetries in the line profiles that hint 
towards velocity fields in a star's atmosphere. Weak emission 
features that need to be clearly distinguished from absorption 
components are promising tools to study circumstellar material. 
Even for observations at lower resolution, the availability of a 
high resolution reference hbrary is required to interpret the re- 
sults. 

We therefore started to build such a library of complete (full 
wavelength coverage from 1 to 5 //m), high S/N NIR spectra 
of representatives of various parts of the HRD using the spec- 
trograph CRIRES (CRyogenic high-resolution InfraRed Echelle 
Spectrograph; Kaufl et al. 2004) at ESO's VLT. With our cho- 
sen setup, it will be possible to study lines fainter than 1% of 
the continuum in strength, and to detect line asymmetries at this 
level. The availability of this library provides an unprecedented 
database for the testing of stellar model atmopheres, the im- 
provement of line Usts of various molecules and atoms, the iden- 
tification of weak lines, and the proper correction of telluric ab- 
sorption features. At this point, an almost complete wavelength 
coverage has been secured for 13 stars. The total library will 
consist of about 25 objects. Further data will become available 
continuously via the dedicated webpage of the project^ The in- 
tention of this paper is to introduce the library to the astronom- 
ical community, discuss aspects of the data reduction, present 
examples, and provide a reference for future exploitation of the 
library spectra. 



2. Observing strategy and instrument settings 

CRIRES is a NIR high-resolution echelle spectrograph operat- 
ing at from 0.96 to 5.3 jum. An adaptive optics (AO) system us- 
ing a curvature wave front sensor can be utilized during obser- 
vations. The instrument works in a near-Littrow configuration. 
Variable entrance slit widths from 0.05 to 3.0 arcsec are offered, 
although in practice only the 0.2 to 0.4 arcsec slits are typically 
used. Active guiding on the slit image in a NIR filter {J, H, or 
K, matching as close as possible the wavelength setting of the 
spectrograph) is part of the standard observing mode. CRIRES 
covers a comparably small wavelength range of a few tens of nm 
per observational setting (depending on the central wavelength). 
To observe the entire range from below 1 to above 5 //m thus 
requires about 200 settings and a considerable amount of ob- 
serving time. A complete scan of a star with K = I mag reach- 
ing a S/N of at least 200 throughout the entire spectral range 
takes almost nine hours and is strongly dominated by observa- 
tional overheads (close to 80%). To achieve our goal of a li- 
brary covering all main spectral types and luminosity classes, 
we limited the hst of possible targets to rather bright sources 
{K < 4?5), and designed the project as a filler programme at the 
Very Large Telescope (VLT) capable of using observing time 
with non-photometric seeing conditions. 

The selection of the targets was driven by the following cri- 
teria: they are (1) representative of a specific spectral type and 
luminosity class or a specific abundance pattern; (2) of sufficient 
brightness in the CRIRES range; (3) of very low rotational veloc- 
ity (typically less than 10 km s"'); (4) found to have no large am- 
pUtude variabihty (less than 1 mag in V); (5) part of the UVES- 



' http://www. univie. ac. at/crirespop 



POP^ sample (Bagnulo et al. 2003) allowing to achieve complete 
wavelength coverage from 300 nm to 5 jum . The latter criterion 
could not be fulfilled in all cases owing to the lack of stars suit- 
able for observations in the NIR range. Furthermore, no suitable 
target exists to represent spectral type O, or Wolf-Rayet stars at 
NIR wavelengths. We included targets with non-standard chem- 
istry, such as a carbon star and a chemically pecuhar Ap star. 
The list of stars observed to date can be found in Table 1, and 
a brief description of each star is given in Sect. 4. Observations 
continue to be performed, and a complete hst of the stars in the 
hbrary will be presented elsewhere at a later date. 

All wavelength settings between 950 nm and 5300 nm of- 
fered by CRIRES were observed, except for those located in- 
between the atmospheric windows, namely 1350-1450 nm, 
1800-1984 nm, 2510-3425 nm, and 4190^620 nm. This pro- 
vided a hst of approximately 200 settings per target. The wave- 
length ranges of the settings overlapped with each other to close 
the gaps between the detector chips. This had the positive side 
effect of increasing the S/N of the final data and allowing us to 
check for instrumental problems. For a few hot stars, we had not 
yet observed the M-band, because we decided to take a closer 
look at the M-band spectrum of one hot star first and identify 
the most promising wavelength settings for observations. This 
would allow us to optimise the use of our observing time. 

For CRIRES-POP, we chose a slit width of 0.2 arcsec, which 
gives a measured resolving power at 2172 nm of 7? = 96 000 
(Smoker 201 1). The AO system was not used because the actual 
spectral resolution is determined by the entrance slit and does 
not depend on the AO correction. 

A pre-dispersion prism and an intermediate slit are used 
within CRIRES to limit the wavelength range to a single order 
before light enters the high-resolution spectrograph, which con- 
sists of a 31.6 hnes/mm, 63.5 degrees blaze echelle grating. The 
science detector is a 4x1 mosaic of Aladdin III InSb arrays. In 
total, 4096x5 12 pixels are read-out; the gap between two succes- 
sive arrays is equivalent to 280 pixels. The detector has several 
defects, particularly a diagonal scratch on the right part of array 

2, a vertical line of bad pixels near the centre of array 4, and a 
larger number of bad pixels at the very left of array 4. Almost all 
of these cosmetic defects will be corrected in the final data set by 
using spectra obtained at different nod positions and overlapping 
wavelength settings. 

There is also an odd-even like pattern present that is aligned 
with the read-out direction. This is in the dispersion direction for 
arrays 1 and 4 and in cross-dispersion direction for arrays 2 and 

3. The effect can be almost completely removed (with residu- 
als in the extracted flux of <k 1%) by flat-fielding and applying 
a correction for detector non-hnearity (see Smoker 2011, Jung 
2010, Uttenthaler 2007). 

At shorter wavelengths (typically below 1.2 /vm), an optical 
ghost is occasionaUy present, which is caused by a reflection 
from the science detector back to the grating. It often does not 
cancel out completely via nodding and is visible as an artificial 
absorption feature in extracted spectra. 

Integration times were chosen in accordance with the 
CRIRES exposure time calculator offered by ESO to achieve a 
S/N of at least 200 in average seeing conditions. Observations 
were always done in an ABBA nodding pattern with typically a 
15" nod throw to account for the background emission. Below 
2.5 fim, total integration times per setting were between 10 s for 



UVES-POP is an electronic library of high-resolution optical spec- 
tra obtained with UVES, the Ultraviolet-Visual Echelle Spectrograph at 
the VLT in the framework of a Paranal Observatory Project. 
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the brightest stars and 120 s for the faintest. In the L-band, total 
integration times were between 60 s and 360 s. For the weaker 
stars, this was split into several integrations at each nodding po- 
sition and up to six nodding cycles. Finally, for the M-band, in- 
tegration times of up to 1200 s were used, again split into sev- 
eral integrations and nodding cycles because of the high thermal 
background contribution. Between 1 and 23 wavelength settings 
were combined into observing blocks of 20 to 60 minutes length. 
We used the no-AO mode. Calibration data were taken during 
day time after each observing night in which our program was 
executed. 

Since CRIRES-POP was designed as a filler programme, the 
various wavelength ranges of a given star were observed over 
up to 1.5 years. A requirement for the selection was, thus, that a 
target was of constant brightness, or only a mildly variable star. 
No airmass limit was set, although we note that only a small 
fraction of observations were obtained at an airmass above 2. 



3. Data reduction 

3.1. First step: standard CRIRES pipeline 

Our CRIRES-POP observations were reduced with recipes avail- 
able from the ESO CRIRES data reduction pipeline 2.0.0 (Jung 
2010). Raw science frames were flat-fielded, and corrected 
for both non-linear detector response and bad pixels. Frames 
from different nodding positions were shifted and combined. 
Finally, the spectrum was extracted and wavelength calibrated. 
The wavelength calibration was determined from either daytime 
ThAr arc lamp exposures for settings below 2.4 //m or sky emis- 
sion lines in the science exposures (above 2.4 fim). We found a 
few instrumental artefacts that in our data and that could not be 
removed during the reduction process. Neighbouring echelle or- 
ders partially overlap at shorter wavelengths, which means that 
part of the spectrum for a given order can be contaminated by ad- 
jacent orders. This is essentially true below 1.3 /urn, where data 
from detectors 1 and 4 cannot be used. When observing spec- 
tral regions with fully saturated telluric features, the actually 
observed line profiles do typically show 1-3 % of residual flux 
in the absorption cores. This is due to in-dispersion straylight 
produced by the CRIRES grating. The eff^ect can be modelled 
quantitatively based on Une-profiles obtained with deep integra- 
tions of neon or krypton arc calibration lamps when the entrance 
slit is illuminated uniformly. When observing a star, i.e. a point- 
like source, the spatial distribution of the scattered light flux in 
the cross-dispersion direction along the sUt then difi'ers from the 
normal spatial profile of the stellar continuum. Thus, the opti- 
mum extraction algorithm employed cannot easily extract pre- 
cise fluxes from the center of absorption lines. Moreover, ow- 
ing to the field-dependent tilt between straylight direction and 
the stellar continuum, there is no easy way to model the ef- 
fect. Hence, the efl'ect cannot be corrected for yet, regardless of 
whether it is caused by telluric absorption or absorption lines in 
the stellar atmosphere. 

Wavelength calibration is a difficult task for CRIRES be- 
cause of the rather small spectral ranges that are measured on 
each detector Wavelength calibration was performed separately 
for each detector and required the availability of a sufficient 
number of calibration lines in the respective spectral region. This 
was not always the case, especiaUy when a ThAr exposure was 
used for calibration. Then a default linear solution was applied. 
The situation, however, might have been improved by using tel- 
luric fines (see below). 



Small differences in the wavelength calibration of the two 
nodded images resulted in a slight degradation of the spectral 
resolution of the combined image that depended on the nodding 
amplitude and specific wavelength setting. We verified on actual 
data that this broadening was smaller than two pixels (i.e. 5.10"^ 
nm/pixel at 1 /zm or 25. 10"^ nm/pixel at 5 yum) for the observing 
strategy used for CRIRES-POP and that the final spectral reso- 
lution was on average degraded by about 10-20% (R=80 000 to 
90 000). Future releases of the CRIRES-POP data will use a new 
CRIRES pipeline recently released by ESO (version 2.1.3), that 
is capable of reducing spectra from A and B nodding positions 
separately. In the future we wiU make use of this feature to re- 
duce spectra where the degradation to the resolution introduced 
by the slight curvature Umits the final data quaUty. 

3.2. Second step: telluric correction using a model of Earth's 
atmosphere 

A general problem in almost all NIR spectra is telluric absorp- 
tion lines. TraditionaUy, they are divided out of a science spec- 
trum using a hot (i.e. feature-poor) "standard" star spectrum. 
In the case of CRIRES-POP, observations of telluric standards 
would have been impractical for two reasons. First, it would have 
led to an enormous amount of caUbration measurements, taking 
into account the large number of wavelength settings involved 
and the need to observe the target and the standard star close in 
time and airmass. Second, telluric standard stars contain a num- 
ber of spectral features, some of them unidentified, such that they 
do not provide ideal 'featureless' spectra and would compromise 
the spectra of our science targets. This is also a reason why we 
included early-type stars in our program, i.e. to explore the in- 
trinsic features of these stars in the IR. Thus, we decided against 
observing telluric standard stars but chose a different approach, 
to model the telluric spectrum at the time of observation. 

Model spectra of the telluric transmission can be calculated 
with freely available radiative transfer codes (e.g. the Line-by- 
Line Radiative Transfer Model [LBLRTM] or the Reference 
Forward Model [REM]), using a database of molecular transition 
parameters (e.g., HITRAN, the High-Resolution TRANsmission 
molecular absorption database) and atmospheric profiles (e.g. 
GDAS^). The modelUng of CRIRES spectra, as described in 
Seifahrt et al. (2010) and Smette et al. (2010), involves the mod- 
elling of the instrumental profile, the scaling of the column den- 
sities of the trace gases in the Earth's atmosphere, as well as a 
wavelength solution for each CRIRES chip as a second or third 
order polynomial. Hence, this approach has the additional ad- 
vantage of providing us with a precise wavelength calibration 
that is often difficult to obtain with standard caUbration tech- 
niques. 

For early-type stars with a clear continuum and few spectral 
features, fitting of teUuric absorption models is semi-automatic. 
Late-type stars with an abundance of stellar lines require a more 
restricted type of modelling and manual intervention caused by 
the strong blending of telluric and stellar features. In extreme 
cases, the instrumental profile and abundance scaling factors 
have to be fixed based on observations of early-type stars in 
the same night and only a solution of the wavelength vector and 
the most variable features (i.e. water vapor) can be obtained for 
the target. Overlapping wavelength regions on different chips 
obtained in neighbouring settings provide an important sanity 
check to the modelling in these cases. 



^ Global Data Assimilation System 

http : //ready . arl . noaa . gov/READYamet . php 
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The current standard of line data of relevant molecules in the 
Earth's atmosphere is set by the HITRAN consortium (Rothman 
et al. 2008). We were able to confront the model spectra with 
our high-quality CRIRES spectra (of hot stars, e.g. t Sco, e Vel), 
to test the quality of the available molecular data. Similar tests 
of molecular data had already been done by some of us for a 
few selected CRIRES wavelength settings (Seifahrt et al. 2010; 
Smette et al. 2010). With the current CRIRES-POP data base, 
however, it becomes feasible to test molecular data for almost 
the entire wavelength range between 0.97 jum and 5.3 fim (a few 
small ranges were omitted because they were affected by heavy 
telluric absorption). We are in contact with the HITRAN consor- 
tium to work out a procedure for a potential improvement of the 
line data. In a first step, the poorly-fit lines will be catalogued. If 
there is a line that deviates from its observed position or strength, 
it can be flagged as suspicious. The HITRAN Ust has quaUty in- 
dicators for each line in the database, that can be used to eval- 
uate whether the observed deviation is expected. If the line is 
consistently off by more than the expected uncertainty listed in 
HITRAN in several observed spectra, this line deserves further 
scrutiny. The wealth of information contained in the CRIRES- 
POP database turns this effort to improve molecular line data 
into a considerable project of its own, which will have a large 
legacy value. Meanwhile, a moderate amount of back-fitting of 
line parameters in terms of position and strength can improve 
the telluric fitting for our project. For the final data set, we aim 
to provide telluric corrections for the entire wavelength range 
covered by our observations where telluric lines are detected at 
the 2 <T level. We expect the estimated performance of the cor- 
rection to be at or above the level shown in Seifahrt et al. (2010), 
especiaUy when the discussed back-fitting of line parameters is 
appUed. 

3.3. Combining individual spectral settings 

Once aU individual one-dimensional (ID) spectra from a single 
star had been extracted (Sec. 3.1) and cleaned of telluric absorp- 
tion features (Sec. 3.2), we combined those spectra into a single 
data product. At this step, we used the spectral overlap of the 
individual settings to improve the final SNR, and to exclude re- 
gions that show contamination from neighbouring orders, optical 
ghosts, or other data defects (e.g. bad pixels). For spectral re- 
gions that include telluric absoption lines, we used the superior 
wavelength calibration delivered by the telluric modelhng and 
combine all spectra after mapping them onto a common wave- 
length vector. At this stage, we also took into account any veloc- 
ity shifts resulting from obtaining the various parts of the spectra 
at different epochs. The combination of several adjacent settings 
into one hour observing blocks ensured that in most cases only 
smaU corrections were needed. For the final Ubrary the spectra 
will be shifted in wavelength to zero heliocentric velocity. 



4. Library content 

In the following sections, we summarised the information from 
the hterature about the global parameters and any determined 
atmospheric parameters for these stars.'' The final release of 
the data will consist of the raw, pipeline-reduced, and telluric- 
absorption-corrected spectra for the complete observed sample. 



The latest version of the CRIRES data for these stars is available 
via http://www.univie.ac.at/crirespop. 



4.1. Current content 

4.1.1. HD 149438 (t Sco), B0.2 V 

At a distance of about 1 50 pc, t Sco (spectral type BO. 2 V) is one 
of the brightest young massive stars of the southern Sco-Cen as- 
sociation. It is therefore one of the most intensely studied early- 
type stars, with observations reported from the X-ray regime to 
the radio. Notable data sets are the pioneering high-resolution 
spectral atlases obtained with both XMM-Newton in X-rays 
(Mewe et al. 2003) and with the Copernicus satelhte in the UV 
(Rogerson & Upson 1977). Unusual stellar wind characteristics 
and unusually hard X-ray emission have been revealed at these 
short wavelengths despite the star's inconspicuousness in the op- 
tical spectrum. This is likely because of pronounced shocks aris- 
ing from its clumped stellar wind (Howk et al. 2000). Moreover, 
the star continues to offer observational surprises, such as the 
detection of its complex magnetic field by Donati et al. (2006). 
tSco has been the subject of numerous quantitative model at- 
mosphere analyses, starting with the first detailed spectral inves- 
tigation of any star besides the Sun, by Unsold (1942). A sum- 
mary of stellar parameters determined in the more recent litera- 
ture is given by Simon-Diaz et al. (2006) and Nieva & Przybilla 
(2008). The most comprehensive non-LTE study to date found 
an effective temperature of T^ff = 32 000±300 K, surface gravity 
log g = 4.30+0.05 (cgs units), microturbulence ^- 5+lkms~\ 
projected equatorial rotational velocity vsini and macroturbu- 
lence ^ of both 4±4kms"', and elemental abundances typi- 
cal of the massive star population in the solar neighbourhood 
(Przybilla et al. 2008), close to solar values. Exceptions are the 
light elements carbon and nitrogen, which indicate a high degree 
of mixing of the surface layers with CN-processed matter from 
the stellar core (Przybilla et al. 2010). 

The first global spectral synthesis modeUng for r Sco of an 
observed spectrum of high quality in the optical was presented 
in Nieva & Przybilla (2011). The observed spectrum was ac- 
quired with FEROS at high resolution (R=48 000) and high S/N 
(~800 in the B-band), and the global spectrum computed as- 
suming non-LTE covers the wavelength region from 3900 Ato 
6850 A. A study of tSco in the NIR shows very good agree- 
ment of parameter determinations in comparison with the analy- 
sis in the optical. The He I A2.05S fim and A2.l \ jum obtained 
with UKIRT/CGS4 (Mountain et al. 1990) could be perfectly 
recovered with non-LTE techniques (Nieva & Przybilla 2007). 
Further quantitative analyses of CRIRES data collected in 2007^ 
and corrected for telluric features following Seifahrt et al. (2010) 
found good agreement between data and models several hydro- 
gen and helium lines, as well as metal lines, when adopting the 
stellar parameters and chemical abundances derived in the opti- 
cal (Nieva et al. 2009, 2011). 

4.1.2. HD 47105 (yGem), AO IV 

Also known as Alhena, y Gem is a nearby AO IV star. Alhena 
has been found to be a wide binary from both spectroscopic 
(Pourbaix et al. 2004) and astrometric (Jancart et al. 2005) mea- 
surements, that has an orbital period of 4614.5 d and an or- 
bital eccentricity of e = 0.89. The Hipparcos measurements (van 
Leeuwen 2007) inferred a distance of 33.5 pc with an uncertainty 
of approximately 2.5 pc. Its rotational velocity was determined 
to be 15 km s"' (Royer et al. 2007), and its heliocentric radial ve- 

The spectra of tSco obtained in the run 079.D-0810(A) have 
been incorporated to the CRIRES-POP database by kind agreement of 
Norbert Przybilla. 



4 



T. Lebzelter et al.: CRIRES-POP 



Table 1. Cross references and spectral types of the stars currently in the CRIRES-POP library, sorted by spectral type. The last column gives the 
year of the observations. 



HD number 


other names 


spectral t5?pe 


RA (2000) 


DE (2000) 


K 


UVES-POP?" 


observing 












(mag) 




periods 


149438 


tSco, HR6165 


B0.2 V 


16:35:52 


-28:12:58 


3.60 


n 


2007"/ 2009-2011 


47105 


7 Gem, HR2421 


AO IV 


06:37:42 


+ 16:23:57 


1.92 


n 


2010 


118022^ 


o Vir, HR 5 1 05 CW Vir 78 Vir 


Alp 


13:34:08 


+03:39:32 


4.88 


y 


2010-2011 


73634 


e Vel,HR3426 


A6II 


08:37:39 


-42:59:21 


3.6 


y 


2009-2010 


20010 


QrFor,HR963, LHS 1515 


F8V 


03:12:05 


-28:59:15 


2.32 


y 


2009-2010 


109379 


jeCrv,HR4786 


G5II 


12:34:23 


-23:23:48 


0.8 


y 


2010-2011 


83240 


1 Leo, HR 3827 


Kl m 


09:37:13 


+06:50:09 


2.66 


y 


2009-2010 


225212 


3Cet, HR9103 


K31 


00:04:30 


-10:30:34 


1.4 


y 


2010-2011 


49331 


HR2508 


Ml 1 


06:47:37 


-08:59:55 


0.6 


y 


2010-2011 


224935 


HR9089, YY Psc 


M3 111 


00:01:58 


-06:00:51 


-0.5 


y 


2009-2010 




Barnard's star, GJ 699, V2500 Oph 


M4 V 


17:57:48 


+04:41:36 


4.52 


n 


2010-2011 


61913 


HR2967, NZGem 


S 


07:42:03 


+ 14:12:31 


0.56 


y 


2010-2011 


134453 


HR5644, XTrA 


c 


15:14:19 


-70:04:46 


-0.6 


y 


2010 



Notes. http://www.sc.eso.org/santiago/uvespop/ Part of the spectrum was obtained within proposal P79.D-0810 and was kindly provided for 
inclusion in the CRIRES-POP library. No M-band spectrum available at the present time. 



locity -12.5 km s-' (Valdes et al. 2004). Zorec et al. (2009) de- 
rived refr=9040+280K for this object using the bolometric flux 
method. Landstreet et al. (2009) used Stromgren and Geneva 
photometry to determine Teff =9150Kandlog^ = 3.5. From high 
resolution optical spectra, they measured v sin ; = 1 1 ± 0.4 km s"^ 
with a microturbulence f = 1.2±0.4kms"'. According to their 
analysis, the iron abundance is close to the solar value. This 
agrees with the results of Adelman & Unsuree (2007), who 
found [Fe/H]=0. 11+0.18. No UVES-POP spectrum is avail- 
able for this star, but it is part of the spectral library pro- 
duced by Valdes et al. (2004) covering the wavelength range be- 
tween 3465 A and 9450 A at a resolution of 1 A FWHM. Flux- 
calibrated low-resolution J, H, and K band spectra of this star 
were presented in the catalogue of Ranade et al. (2007). 

4.1.3. HD 118022 (o Vir), Alp 

o Vir (Alp) has the distinction of being the star for which 
the pioneering work of Babcock (1947) employed the Zeeman 
eflBCt to determine the magnetic field. Its magnetic field has 
since been well-studied (mean magnetic field < > = 3.0 kG, 
Ryabchickova et al. 2008), in addition to other atmospheric phe- 
nomena. An effective temperature of 9460 K was determined by 
Netopil et al. (2008) in good agreement with the value used in 
the analysis of Ryabchickova et al. (2008). The latter also found 
log ^ = 4.0. Its light curve displays a small amplitude variation 
(0.07 mag in V), and it is classified as an CVn type vari- 
able, its variability being caused by the rotation of an inhomo- 
geneous stellar surface. Its measured value of v sin j is 10 km s"' 
(Ryabchikova et al. 2008), which we verified by performing a 
synthetic spectrum analysis of UVES and CRIRES data. The 
optical spectrum displays a richness caused by the rare earth el- 
ements, which appears clearly evident at NIR wavelengths as 
broad but shallow features. The presence of Eu ii near 1 jxm was 
noted by Wahlgren (2011). The IR spectrum will be utilised to 
identify atomic fines from heavy elements, and lines for mag- 
netic field studies in chemically pecuUar stars. 

4.1.4. HD 73634 (eVel),A6ll 

The bright giant e Vel (HR 3426) was proposed as a low v sin ; 
A6 11 spectral standard by Gray & Garrison (1989). The UVES- 



POP data illustrate the sharp-lined nature of the spectrum, along 
with several broader, weak emission features at red wavelengths. 
From our synthetic spectrum fitting of the UVES-POP spec- 
trum, V sin i was found to be 8.5±0.2 km s"' . Basic data derived 
from standard calibrations were presented by Snow et al. (1994), 
namely T^s = 8100 K, Mboi = -2.8, and rv = 19 kms'^ 

4.1.5. HD 20010 (a For), F8 V 

Saffe et al. (2008) derived values for Teff and \ogg of 
6094 K (6226 K) and 4.09 (3.77) dex, respectively (depending 
on the Teff-log ^-colour relation used). According to these val- 
ues, the metaUicity [Fe/H] is -0.64 (-0.62). A more compre- 
hensive spectroscopic analysis by Gonzalez et al. (2010) found 
reff = 6170±35K, logg = 3.80±0.08 and [Fe/H] = -0.21 ±0.03, 
and Bruntt et al. (2010) derived Jeff = 6015 K, logg = 3.80, 
vsinj = 3.9kms-\ ^ = 3.7kms-\ and [Fe/H] = -0.28, with flie 
latter value also being representative of many more chemical 
species. A parallax of 70.24 mas was derived from Hipparcos 
observations (van Leeuwen 2007), suggesting a distance of only 
14pc. Oudmaijer et al. (1992) searched for possible dust around 
this star using IRAS data, but found only a very small IR excess. 

4.1.6. HD 109379 (/3Crv), G5 II 

The most recent determination of the parameters of this G5 

II star was presented by Lyubmikov et al. (2010), giving 
Teff = 5100±80K, log^ = 2.52±0.03, and [Fe/H] = +0.15. Using 
the Hipparcos parallax (van Leeuwen 2007), a distance of 44.6 
pc was derived, and amass of 3.7±0.1 Mq determined. The study 
of Takeda et al. (2008) found similar values of Teg- = 5145 K, 
logg = 2.56, and [Fe/H] = -0.01, in addition to log L/Lq = 2.23 
and M=3.31M0. Older investigations (e.g. Luck & Wepfer 
1995, Brown et al. 1989) came to similar conclusions. The he- 
liocentric radial velocity was measured to be -8kms"' (e.g. 
Buscombe & Morris 1958), and vsmi was around 4kms"' 
(Gray 198 1, Fekel 1997). A radius of 1 1 /?© is listed in Pasinetti 
Fracassini et al. (2001). 

While the C/O ratio was found to be close to the solar value, 
the two elements C and O seem to be depleted in this star relative 
to the solar value according to the analysis of Luck & Wepfer 
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(1995). An earlier claim of an enrichment in barium (Keenan & 
Pitts 1980) could not be confirmed (Eggen & Iben 1991). 



4.1.7. HD 83240 (10 Leo), K1 III 

Stellar parameters for this Kl giant, also known as the bright 
star HR3827 or 10 Leo, were derived by Soubiran et al. (2008). 
They found Tes and log^ to be 4682 K and 2.45 dex, respec- 
tively. Soubiran et al. evaluated a metallicity [Fe/H] of -0.02. 
Compared to Arcturus, which is thought to be a metal-poor ob- 
ject (e.g. Peterson et al. 1993), we thus have here a K giant ref- 
erence spectrum of solar metalUcity. The Hipparcos distance is 
69 pc. The star is a spectroscopic binary with an orbital period 
of 2834+4 d and y=+20.0±0. 1 km s"' (Griflin 1985). No trace of 
the secondary star was detected, the orbital parameters indicate 
that it is of low mass. Pasinetti Fracassini et al. (2001) listed a 
radius of 14Rq for this object. No indication of light variabiUty 
was found (Percy 1993). 

Mishenina et al. (2006) measured abundances of vari- 
ous key elements from the spectrum of HD 83240. They 
found that log A(C) = 8.25, log A(N) = 8.2, log A(0) = 8.75, and 
logA(Li) = 1.1. This means a slight underabundance in carbon 
and an overabundance in nitrogen, which agrees with a post first- 
dredge-up red giant. Silicon, calcium, and nickel are solar-like in 
abundance. In a second paper, Mishenina et al. (2007) provided 
abundances for various neutron-capture elements. None of the 
element abundances differed from the solar value by more than 
the typical error given by the authors. 



4.1.8. HD 225212, K3 I 

The most recent and most detailed study of the global parame- 
ters of this K-type supergiant was provided by Smiljanic et al. 
(2006), finding Teff = 4052K, log § = 0.75, ^ = 2.95 km s^^, and 
[Fe/H] = +0. 1 . They derived Mboi = -4.42, and estimated a mass 
of this object of 7.0-7.9 M©. Their heliocentric radial velocity 
value of -42kms~' is in excellent agreement with other mea- 
surements in the literature (De Medeiros et al. 2002, Valdes et 
al. 2004, Eaton & Williamson 2007). Earlier studies of the stel- 
lar parameters found quite diverse results: van Paradijs (1973) 
measured a similar Teff =4100 K, but a general overabundance 
of heavy elements by a factor of 1 .6, while Luck & Bond (1980) 
derived for almost the same effective temperature a metallicity of 
[Fe/H] = -0.2. This scatter may result from the different values 
of the microturbulence used. Lower values for the temperature 
were suggested by Blum et al. (2003, 3860 K) and MalUk (1999, 
3689 K). Values very similar to the results by Smiljanic et al. 
(2006) can be found in the Ubrary of Sanchez-Blazquez et al. 
(2006). 

Smiljanic et al. derived abundances of [C/Fe] = -0. 14, 
[N/Fe] = -1-0.41, and [O/Fe] = -0.02. These findings agree with 
predictions for an atmospheric composition past first dredge- 
up for a non-rotating star. Consistently, a rotational velocity of 
5.8±1 km s"^ has been measured (De Medeiros et al. 2002). 

No indications of binarity have been detected (Eggleton & 
Tokovinin 2008). Koen & Eyer (2002) found evidence for a two 
hour variability with an amplitude of a few imnag from the epoch 
photometry of the Hipparcos catalogue. The parallax of this ob- 
ject is 1.56+0.31 mas (van Leeuwen 2007). 



4.1.9. HD 49331, Ml I 

Two studies were dedicated to the measurement of stellar pa- 
rameters of this M supergiant. Smith & Lambert (1986) de- 
rived Teff = 3600 K, logg = 0.7, and Mboi = -3.9. The Hipparcos 
distance leads to a somewhat brighter Mboi - -4.5 (using the 
V brightness and the bolometric correction of Worthey & Lee 
2011). Smith & Lambert also determined the abundances of 
various elements and found [Fe/H] = +0.14, '^C/'^C = 18, and 
^^C/^^O = 0.32, with no indications for an enhancement in s- 
process elements. The more recent investigation of Wylie-De 
Boer & Cottrell (2009) confirmed temperature, log g value, and 
the slightly super-solar metallicity of HD 49331, but measured 
a higher C/O ratio of 0.54 and a mild enhancement of s-process 
elements. The latter result is in agreement with an early spectral 
classification of the star as MS due to enhanced bands of YO 
and ZrO (Yamashita 1967). However, clear signs of technetium, 
a key indicator of recent third dredge-up, are missing (Little et 
al. 1987). 

A mean radial velocity of 26.7 km s ' was measured by 
Brown et al. (1990). 

4.1 .1 0. HD 224935 (YY Psc), M3 III 

The bright, cool giant YY Psc (M3 111) has an accurate paral- 
lax determined by the Hipparcos satelUte of 7.55±0.59mas (dis- 
tance of 132pc, van Leeuwen 2007). According to the GCVS 
(Samus et al. 2007-2010), it displays irregular variability with 
an amplitude of 0.1 mag at visual wavelengths, which was con- 
firmed by Hipparcos. Owing to the small ampUtude of variabil- 
ity, we assumed that variability had no impact on the analy- 
sis of the CRIRES spectrum. The angular diameter of YY Psc 
was measured both via long baseline interferometry (Dyck et al. 
1998) and by the lunar occultation technique (Fors et al. 2004). 
Both methods provided a value close to 7 mas. Dyck et al. also 
derived Teff = 3647+184 K, and a similar temperature was found 
in the investigation of Dumm & Schild (1998), who listed a ra- 
dius of 109 Rq and a mass of 2.5 M©. Studies by Sloan & Price 
(1998) and Kwok et al. (1997) found YY Psc to be dust-free. 

Eaton & WilUamson (2007) measured a radial velocity of 
12.05+0.43 km s~' from a series of 103 observations. The rota- 
tion velocity, vsin;, was determined by Zamanov et al. (2008) 
to be either 2.5 or 3.9 km s~\ depending on the method ap- 
pUed. Lazaro et al. (1991) measured a '^C/^^C ratio of ten from 
low dispersion IR spectra, with a carbon abundance that is defi- 
cient by [C/HJ = -1.25, relative to the Sun. Using the CRIRES- 
POP spectrum for a preliminary synthetic spectrum analysis with 
a model of the parameters Jeff = 3500 K and logg = 1.0 (cgs), 
Wahlgren et al. (2011) confirmed the >^C/'^C value of Lazaro, 
but not the carbon abundance, which was found to be near solar 
in value. The abundances of both nitrogen and oxygen (both near 
solar) were also determined and enhancements for the heavy el- 
ements (Co, Ge, Sr, Cs, Ba, and Eu). Lines of Si, Y, and Zr were 
also identified in the CRIRES spectrum. The abundance of iron 
is determined to be solar. The value for v sin i is 4 km s"' . 

4.1 .1 1 . Barnard's star (V2500 Oph, GJ 699), M4 V 

Probably the most prominent star in the CRIRES-POP library 
is Barnard's star (Barnard 1906), the star with the largest 
known proper motion. It is the brightest known M dwarf 
and the only observable M dwarf within our CRIRES-POP 
brightness limits. Stellar parameters found in the litera- 
ture for this star agree very well. The most recent study of 
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Chavez & Lambert (2009) measured reff=3134K, log g = 5.1, 
^ = 0.6 kms-\ [Fe/H] = -0.8, and [Ti/H] = -0.61. Using 
flux ratios from various wavelength bands, Casagrande 
et al. (2008) found Teff = 3145±69K, and Segranasan et 
al. (2003) derived Teff = 3163±65K, log g = 5.05+0.09, 
M = 0.158+0.008 Mq, and = 0.196+0.008 /^o using VLTl 
measurements. Dawson & De Robertis (2004) determined a 
luminosity of L = (3.46+0.17)10"^Lo and Teff = 3134+102 K. 

The star's radial velocity was found to be -110.5 km s"' 
(Nidever et al. 2002) with a very high degree of stabihty at the 
3-5 ms~^ level (e.g. Bean et al. 2010). Owing to its proximity 
and large space motion, it shows a non-negligible amount of 
secular acceleration of approximately 4.5 m s"' yr"' (Kuerster et 
al. 2003). Benedict et al. (1998) measured a rotation period of 
approximately 130 days for this star. The rotational velocity is 
thus very small and only upper limits have been reported (e.g., 
V sin ;• «:2.8 km s"^ , Mohanty & Basri 2003). This old object (7- 
1 3 Gyr, Feltzing & Bensby 2008) shows only a low level of mag- 
netic activity (e.g. Hiibsch et al. 1999), although Paulson et al. 
(2006) investigated a rare flare event on this star. 



4.1.12. HD 61913, (NZ Gem), S 

The spectral classification of this star is slightly unclear It 
was classified as M3S by Keenan (1954). Yamashita (1967) 
found that ZrO and YO were slightly enhanced, while Smith 
& Lambert (1988) found no enhancement in s-process ele- 
ments. Nevertheless, the star is listed in several later papers 
among MS/S stars, and our own inspection of the UVES- 
POP spectrum indicates a spectral type of MS owing to the 
clearly visible ZrO bands. Famaey et al. (2005) sunamarized 
data about this star's distance and kinematics, namely d=324pc 
and rv/,e;=-15.2kms"'. This velocity is in good agreement with 
findings from Jorissen et al. (1998). The latter authors also 
found some velocity jitter in this star, but no clear indication 
for binarity. Van Eck et al. (1998) derived an absolute bolo- 
metric magnitude of -4?01. Cenarro et al. (2007) calculated 
Teff =3530 K and log g = 0.7 from the spectral type and lumi- 
nosity class (M3II-III). NZ Gem shows no lines of technetium 
(Lebzelter & Hron 2003), which suggests that any possible en- 
hancement in s-process elements is due to an extrinsic enrich- 
ment. However, no indications for a binarity in this object were 
found in the study of Eggleton & Tokovinin (2008). 
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Fig. 1. Location of the CRIRES-POP stars in a HRD. Parameters taken 
from the literature or derived as described in the text. Filled symbols 
mark stars presented here, and open symbols the forthcoming objects. 
The solid lines indicate evolutionary tracks for solar-metallicity models 
of 10, 5, 3, 2, 1, and 0.6 Mq taken from Bertelli et al. (2008, 2009). 



spectively. According to Kipper's analysis, X TrA seems to be 
slightly metal-poor ([Fe/H]= -0.5) with a supersolar abundance 
of the s-process elements Y and Zr. The carbon isotopic ratio 
'^C/'^C was found to be 31. A C/O ratio of l.I was used for all 
carbon stars investigated by Kipper. The star probably displays 
lines of Tc (Little et al. 1987), which we confirmed by inspecting 
the UVES-POP spectrum. Together with enhanced abundances 
of Y and Zr, this suggests that the star is an intrinsic AGB star, 
although the carbon isotopic ratio reported in the literature seems 
a bit low. Winters et al. (2003) derived a mass loss rate of about 
5x10"^ Moyr"^ from thermal CO lines. 

4.2. Forthcoming targets 

The CRIRES-POP observations are expected to continue for an 
additional year to reach a total sample size of approximately 25- 
30 targets. Spectra have been partly obtained for seven further 
targets, which are listed in Table 2. Information about the stars 
will be published onUne* as soon as observations are completed. 



4.1.13. HD 134453 (X TrA), C 

The group of C-rich stars is represented in our library by the 
bright star HR 5644 (X TrA). A spectral type C5,5 was assigned 
to this object (Wallerstein & Knapp 1998). The star is an irreg- 
ular variable with a maximum amplitude of 1 magnitude in the 
visual according to the GCVS (Samus et al. 2007-2010). Percy 
et al. (2009) estimated a variabihty period of 500 days. In a mul- 
tiplicity study of bright stars by Eggleton & Tokovinin (2008), 
the star is listed as a single object. A distance of 460 pc, based on 
the Hipparcos parallax, is given by Winters et al. (2003). These 
authors also list a LSR- velocity for X TrA of -2.5 kms~\ cor- 
responding to a heUocentric velocity of +2.5kms~'. This is in 
agreement with measurements by Schoier & Olofsson (2001). 

Kipper (2004) summarized earlier estimates of the star's 
effective temperature and absolute bolometric magnitude, 
which he used to derive surface gravity and microturbulence 
from the UVES-POP spectrum. He determined the values of 
Jeff = 2700 K, Mboi = -6.0, log^ = -0.8, and^ = 2.3kms ^ re- 



4.3. The CRIRES-POP HRD 

We include in Fig. 1 the location of the stars of the CRIRES-POP 
library in a Hertzsprung-Russell diagram (HRD). Temperatures 
and luminosities were taken from the literature, and for each 
source we refer to the individual target descriptions in the para- 
graphs above. If no luminosity had previously been calculated, 
we used the Hipparcos parallax (van Leeuwen 2007), tem- 
perature, colour, and V magnitude from the literature as de- 
scribed above, and the bolometric correction from Worthey & 
Lee (2011) to estimate a value for Mboi. Luminosities from the 
literature were converted into Mboi using Mboi,© = 4?74. In Fig. 1 
we did not included the star HD 73739 owing to a lack of suffi- 
cient information about this object in the literature. 

We include in Fig. 1 evolutionary tracks for stars with solar 
metallicity and masses of 10, 5, 3, 2, 1, and 0.6 M©, taken from 
Bertelli et al. (2008, 2009). The tracks are primarily for illustra- 
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Table 2. Forthcoming additions to the CRIRES-POP library, sorted by spectral type. 



HD number other names 


spectral type 


RA (2000) 


DE (2000) 


K 

(mag) 


UVES-POP? 


120709 


3 Cen A, HR 5210 


B5 IIIp 


13:51:50 


-32:50:39 


4.97 


y 


80404 


L Car, HR 3699 


A8I 


09:17:05 


-59:16:31 


1.53 


y 


146233 


18 See, HR 6060 


G2 V 


16:15:37 


-08:22:10 


4.19 


n 


99648 


T Leo, HR 4418 


G8I 


11:27:56 


+02:51:23 


2.83 


y 


138716 


37 Lib, HR 5777 


Kl V 


15:34:11 


-10:03:52 


2.24 


y 


209100 


eInd,HR8387 


K5 V 


22:03:22 


-56:47:09 


2.20 


y 


73739 


MN Vel 


M7II/in 


08:38:01 


-46:54:15 


-0.19 


n 



live purposes. Owing to the inhomogeneous collection of tem- 
perature and luminosity values presented here, we did not aim 
to perform fits of the mass and evolutionary status of the vari- 
ous targets. Our intention was instead to highlight the parameter 
range covered by CRIRES-POP. 

5. Results and discussion 

5.1. A first glance at the catalogue 

As outlined above, the final CRIRES-POP data products should 
clearly have an accuracy that exceeds the average level of the 
standard CRIRES pipeline products. In particular, the wave- 
length solution has been improved by the usage of telluric lines 
(see Sect. 3.2), cosmetic defects have been corrected where pos- 
sible (see Sect. 2), and the full spectral resolution with a re- 
solving power of R ^ 100000 will be maintained, compared 
to /? » 80 000-90000 from the standard pipeline reduction (see 
Sect. 3.1). We started the production of the final library, and pre- 
sented here the first examples. In Fig. 2, we illustrate the process 
of telluric correction and fine tuning of the wavelength cahbra- 
tion briefly described in Sect. 3.2. This method has been applied 
in three representative wavelength regions for the 13 stars listed 
in Table 1 . We have selected these parts of the spectra because 
they include the positions of some well-known hydrogen lines, 
several atomic fines, and a large collection of molecular fines, 
i.e., interesting features that can be seen in all stars of our Ubrary 
in this presentation. 

Figures 3 to 15 show the spectra of the three wavelength 
regions for each star in the library and provide line identifica- 
tions. In our presentation, we follow the sorting given in Table 1 . 
Heliocentric velocity corrections have been appUed to all spec- 
tra. The accuracy of the wavelength solution, derived with the 
help of the positions of the telluric lines, is better than one third 
of a pixel, and for the hot stars the accuracy is close to one tenth 
of a pixel. A representative synthetic telluric spectrum for these 
regions can be found in Fig. 16. The telluric corrections appUed 
to the selected spectra shown here leave residuals in the cores of 
the most prominent telluric lines of ~1% of the continuum flux 
for all but one noticeable exception at 1278.5 nm. The correc- 
tions were achieved without back-fitting of line parameters and 
using a single Gaussian instrumental profile. The figures clearly 
illustrate the high quality of the CRIRES-POP data. 




1.0 
0.9 
0.8 
0.7 
0.6 
0.5 



pipeline output_ 
I I 1 



r]if\f\f\ 



T 




telluric model 
I I 1 




final spectrum _ 
J 1 I 1 



4040 4045 4050 4055 4060 

wavelength [nm] 



4065 



Fig. 2. Example of our method of correcting for the telluric absorption. 
The top panel shows the spectrum of HD 73634 (A6 11) around the Br-or 
line at the level of the CRIRES pipeUne data product after merging of 
two slightly overlapping spectral regions. The middle panel shows the 
telluric model calculated to fit the telluric lines in the observed spec- 
trum. The bottom panel shows the resulting spectrum after division by 
the telluric spectrum and the application of some basic cosmetic correc- 
tions. Note that the emission core in the broad hydrogen line is real and 
not a residual of a telluric line. 



5.2. Line identifications 

In the presented spectra, atomic and molecular features have 
been identified. We used the identifications given in the infrared 
atlas of Arcturus (Hinkle et al. 1995), the fine lists provided by 
Kurucz & BeU (1995), the fist of iron lines in Nave et al. (1994), 
the HITRAN database from Rothman et al. (2008), and the iden- 
tifications of SiO and OH fines at 4 fim given in Ridgway et al. 



(1984). The attribution of atomic or molecular transitions to the 
observed features was based on identifications for similar stars 
presented in the literature or the most probable origin of the line, 
and have in some cases been verified by spectral modelling. One 
of the weakest lines identified is probably a Na line at 4044.5 nm 
in HD 73634 (Fig. 6), which has a central depth of not more than 
2 %. Unidentified lines even weaker than that can be seen e.g. in 
the //-band spectrum of the same star. 
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Line identifications for t Sco were taken from Nieva et al. (in 
preparation). Non-LTE population inversions in the outer atmo- 
sphere give rise to emission cores in Br a in the hotter stars (see, 
e.g., Kaufl 1993, Przybilla et al. 2004), that are most notable for 
tSco, Pa/3 also displaying the phenomenon. 

Line identification in the spectra of the majority of the 
observed targets was straightforwardly achieved; however, the 
identifications for two stars (the M dwarf, Barnard's star, and 
the carbon star X TrA) required a special treatment owing to the 
lack of information available in the literature for these spectral 
types. For these two stars, the spectra and line identifications are 
discussed in detail below. 



5.2.1 . Barnard's star (GJ 699) 

The closest high-resolution template to our spectrum of 

Barnard's star (GJ699), is the FTS atlas spectrum of a sunspot 
(i.e. dark umbra) acquired by Wallace et al. (2001, 2002). 
However, the available spectra of the sunspot and GJ 699 differ 
remarkably, which is not completely surprising, given the high 
magnetic field and hence low gas pressure in a sunspot (e.g., 
Solanki 2003), mimicking a low surface gravity environment 
(Amado et al. 2000). In addition, GJ 699 has a lower metallic- 
ity ([Fe/H] = -0.5, Jones et al. 2002) compared with the Sun. 

Moreover, the high magnetic field in the sunspot leads to 
strong Zeeman splitting and broadening of atomic (e.g. Ti I and 
Ca I around 1250 nm) and molecular (FeH) lines, which we did 
not detect inGJ699. 

The FeH molecule is a strong opacity source in the red and 
NIR for low mass stars, and those lines were readily detectable in 
our J and H band spectra of GJ 699. Lines were identified based 
on fine lists published in Dulick et al. (2003) and Hargreaves et 
al. (2010), where only the strongest lines are marked in the plots 
shown in Fig. 13. FeH is magnetically sensitive and suspected to 
weaken with lower surface gravity. The Zeeman broadening and 
the low-gravity environment thus greatly suppress those lines in 
the sunspot spectra, where they are mostly absent in the spectral 
region analysed here. To the contrary, low surface gravity en- 
hances the line depth of CO, and we found CO second overtone 
fines in the sunspot spectra around 1650 nm, which were absent 
inGJ699. 

The 4050 nm spectrum of GJ 699 is dominated by a forest 
of fines, which we identified as water vapor. In contrast to FeH, 
where we identified all lines directly with their counterpart in the 
line fists, we could only provide a tentative identification of the 
H2O lines because the latest line lists of Partridge & Schwenke 
(1997) and Barber et al. (2006) provided only a rough match. 
A number of lines in our spectrum seem to be Umited by the 
predicted line positions for the same transition in both fists and a 
number of predicted lines appear to be missing, while other lines 
found in the GJ 699 spectrum have no unambiguous counterpart 
in those lists. All H2O lines found in GJ 699 are present in the 
FTS sunspot spectra (although much weaker than in GJ699), 
whereas they remained mostly unidentified in the Wallace et 
al. atlas, except for a few lines that are indeed identified with 
H2O. We note that the main opacity source in the sunspot spec- 
tra around 4050 nm are SiO lines (some isolated 2-0 lines and the 
3-1 bandhead, see, e.g., the spectrum of NZ Gem, Fig. 14), that 
are, however, extremely weak in GJ699. The same is true for a 
small number of OH (4-1) lines falling in this window. The Bra 
fine, present in the sunspot spectra, was undetected in GJ 699. 



5.2.2. HD 134453 (XTrA) 

The second star where we could not rely on identifications of 
even the stronger lines in the fiterature was the carbon star 
HD 134453. High resolution NIR studies on C-stars have lead 
to interesting results on, e.g., the abundance of fluorine (Abia et 
al. 2010, Jorissen et al. 1992), but a systematic investigation of 
the spectral content is stiU missing. 

In the / and H bands, lines of CN and C2 dominate the spec- 
trum. Blending is a severe problem, hence isolated lines are rare. 
To determine the contribution of the various molecules to an in- 
dividual spectral feature, we calculated synthetic spectra using 
the COMARCS models of Aringer et al. (2009). The fine lists of 
CN and C2 that we used were those from J0rgensen (1997) and 
Querci et al. (1974), respectively. For the model parameters we 
applied the values of Kipper (2004). Unfortunately, the molec- 
ular data in the fine fists lack sufficient precision for this high 
resolution study, hence the synthetic spectrum differs strongly 
from the observations (see also the discussion in Lebzelter et al. 
2008 or Lederer et al. 2009). 

Focusing on the observed features, which were at least 
roughly reproduced by the synthetic spectrum, we re-calculated 
the model spectra by removing one molecule (CN, C2, CO) at a 
time. In this way, we were able to attribute a number of lines or 
parts of line blends to a specific molecule. However, a significant 
fraction of features could not be identified unambiguously. For 
the C2 lines (Phillips system), the line list did not include identi- 
fications of the transitions, thus it is not provided in Fig. 15. 

The model was similarly unable to reproduce the section of 
the spectrum around 4 fim. Wallace & Hinkle (2002) identified 
bands of CH and CS in medium-resolution spectra of carbon 
stars in this wavelength range. We were able to attribute series 
of lines in our spectrum to transitions of the 2-0, 3-1, and 6-4 
bands of '^CS. For this analysis, we used the line lists of Winkel 
et al. (1984) and Chandra et al. (1995). Two likely identifications 
of atomic lines are also given. 

5.3. Impact on atomic and molecular line lists 

Even at this preliminary stage, the potential of the CRIRES- 
POP fibrary is clearly evident for an accurate testing of existing 
atomic and molecular line fists. Laboratory and theoretical inves- 
tigations to create atomic data have in the past focused on the op- 
tical and ultraviolet wavelength regions to support ground-based 
and space spectroscopy. Nevertheless, the infrared wavelength 
region between approximately ljum and 5 /im contains numerous 
atomic transitions. Databases lack transition data for many ions 
at IR wavelengths, giving the false impression that this region 
is uninteresting for studies that are compelled to utilise atomic 
lines. 

Prior to the inauguration of the CRIRES and the Phoenix 
(Hinkle et al. 1998) instrument, the low spectral resolution of 
earlier instruments did not necessitate accurate wavelengths or 
transition probabifities for atomic fines owing to the effects of 
fine blending^. In the presence of molecular features typical of 
cool stars, atomic lines tend to become incorporated into the 
broader features. The incomplete treatment of terrestrial atmo- 
spheric fines contributes a similar eff'ect. 



' A limited amount of high spectral resolution data from instruments 
such as the Kitt Peak National Observatory Fourier transform spectro- 
meter and the solar spectrum obtained from ATMOS (Geller 1992) and 
the ACE satellite (Hase et al. 2010) are available and we clearly need to 
produce accurate atomic and molecular data. 
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The vast majority of data for atomic transitions at IR wave- 
lengths originates from atomic structure calculations. These 
studies were incomplete in terms of the number of transitions 
and possess accuracies that remain, in many cases, undeter- 
mined. Modern experimental techniques for determining oscil- 
lator strengths from transition probabilities have been based on 
the measurement of energy level lifetimes and branching frac- 
tions from line intensities (Wahlgren 2010, 201 1), but few NIR 
Unes have been analysed for atomic parameters. 

Fortunately, many IR transitions originate from energy levels 
that are studied experimentally for optical and ultraviolet transi- 
tions. Additional lifetime measurements will no doubt need to be 
made to enable the strongest and/or best placed lines in the NIR 
of various elements and ions to be analysed. The large amount 
of ongoing laboratory analyses should therefore be targeted to- 
ward recording emission Une spectra for intensity measurements 
to determine transition branching fractions. The wavelength re- 
gion between 0.9 yum and 2.5 yum is of particular importance for 
atomic data, as this region contains known transitions of many 
post iron-group elements. The iron-group elements are of inter- 
est at all wavelengths to provide the wavelength standards that 
can be used to adjust the CRIRES wavelength scale in regions 
where few calibration lines exist, and to apply accurate oscilla- 
tor strengths to place abundance analyses at UV, optical, and IR 
wavelengths on a similar scale. Finally, the spectral library wiU 
provide valuable comparison data for testing line-broadening 
theories. Few precise data are available so far, but they are in 
high demand for quantitative analyses using not only the hydro- 
gen and helium lines, but also several rather strong metal transi- 
tions. 

6. Conclusions and outlook 

We have presented the first results and relevant background in- 
formation for our high-resolution spectral library in the infrared 
called CRIRES-POR The spectral data that we have presented 
illustrate the potential of our library not only for the studies of 
individual stars that are underway within the CRIRES-POP con- 
sortium, but also comparative studies spanning a range of spec- 
tral types. As an immediate highlight, we considered the case 
of Barnard's star. Our high-resolution spectrum has allowed for 
the first time a direct identification of photospheric water vapor 
lines in a cool M dwarf and demonstrated the incompleteness 
of present-day H2O line-lists. The large wavelength coverage of 
the CRIRES-POP observations, combined with its high spectral 
resolution, ensures that they provide an invaluable resource for 
prioritising laboratory measurements of both wavelengths and 
transition probabilities. The diversity in the properties of stars 
observed within the CRIRES-POP program will provide spec- 
tral lines with a large range of excitation energies and ionisation 
stages. We invite and encourage the scientific community to use 
the reduced data of our Ubrary for further studies of the individ- 
ual stars, but also for comparison, model testing, and preparation 
of observations, but kindly ask researchers to give reference to 
this pubUcation. 
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Fig. 3. Three representative sections of tiie CRIRES-POP data of HD 149438 (B0.2 V). Telluric lines have been removed as described in the text. Known atomic and molecular features are marked. 
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Fig. 4. Same as Fig. 3 for HD 47105 (AO IV). 
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Fig. 5. Same as Fig. 3 for HD 118022 (Alp). 
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Fig. 6. Same as Fig. 3 for HD 73634 (A6 II). 
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Fig. 7. Same as Fig. 3 for HD 20010 (F8 V). 
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Fig. 16. Example of a synthetic telluric spectrum used for the correction of CRIRES-POP data. 



